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ABSTRACT 


Maiomeeee, Grag, and inertia coefficients have been 
determined experimentally for various tube bundles and for 
two cylinders at various relative positions subjected to 
harmonically oscillating flow. The force coefficients for 
the inline force have been analyzed through the use of 
Morison's Equation and Fourier Analysis. The transverse 
force or the lift force has been expressed in terms of its 
maximum value. 

The results have shown that the interference between the 
cylinders can give rise to complex flow patterns and to un- 
expectedly large force transfer coefficients. The results 
with two-cylinder experiments have shed considerable light 


on the wake interference in harmonic flows. 
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NOMENCLATURE 


Drag coefficient 
Inertia coefficient 
Diameter of individual cylinder comprising riser 
» Ds 
1. 
fey)” 
Length of the cylinder 
Total force exerted on the cylinder by the fluid 
Measured force 
Instantaneous velocity, U=-Um cos wt 
Maximum velocity in a cycle 


Frequency of vortex shedding 


Period of flow oscillation 











eg, tk 
Keulegan-Carpenter Number, K = - 
Oe D 
Reynolds Number, Re = . 
r., D 
Strouhal number, St = ij 
m 


Lift coefficient 

Pere LOrce 

Period of flow oscillation 

Time 

Angular frequency (a2) 
2 Re 

Frequency parameter, 8 = D /vT = z= 

einemarle viscosity Of £luid 


Fluid density 


Angle between two cylinder-center line and flow 
alsect: Lon 
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Na eODUCTION 


A body's resistance to flow is strongly affected by what 
Surrounds it. When two bodies are in close proximity, not 
only the flow about the downstream body but also that about 
the upstream body may be influenced. Examples include con- 
denser tubes in heat transfer; columns in pressure suppression 
pools of nuclear reactors; risers, piles, and other tubular 
Structures in off-shore engineering; turbine and compressor 
blades in mechanical or aerospace engineering; and high-rise 
Peaidings, cooling towers, and transmission lines in civil 
engineering. The quantification of the interference effects 
in terms of the pressure distribution, lift and drag forces 
On individual members, vortex shedding frequency, and the 
dynamic response of members of the array in terms of the 
governing flow and structural parameters constitute the 
essence of the problem. 

There are infinite numbers of possible arrangements of 
two or more bodies positioned at right or oblique angles to 
the approaching flow direction. In wavy or time-dependent 
flows one needs the lift, drag and inertia coefficients for 
all members of the array. Evidently, the members of the 
array may not be all parallel and normal to the flow. The 
quantification of the flow interference on lift, drag and 
inertia coefficients for cylinders with relative inclinations 


and spacings in a design-wave environment is an exceedingly 
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complex problem. In the absence of data on the interaction 
between drag and inertia coefficients in wavy and harmonic 
flows for cylinder groups one is tempted to use steady-flow 
results for the drag coefficient and the unseparated potential 
flow results for the inertia coefficient. The inertia coef- 
ficient for a group of cylinders in inviscid unseparated 

flow may be obtained through the use of the method of images 
yl, 2 Th or through the use of the linear potential theory 
including wave diffraction Ties, 4 /. These analyses do not 
deal with the effects of separation and vortex shedding. Con- 
sequently, the results are more appropriate to the determina- 
tion of earthquake forces and wave forces on large bodies 
rather than to the evaluation of the inertial components of 
the force in the drag/inertia dominated regime. 

It is evident from the foregoing that experiments must be 
Carried out with carefully selected tube bundles and cylin- 
der combinations in order to develop some understanding of 
the flow interference and to provide data for body combina- 


tions of practical importance. 
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IIT. REVIEW OF PREVIOUS INVESTIGATIONS 


A careful review of flow interference between two circular 
cylinders in various arrangements in steady flow has been pre- 
sented by Zdravkovich / 5 / where an extensive list of re- 
ferences may be found. Numerous studies have shown that the 
changes in drag, lift and vortex shedding are not necessarily 
Pem@emiuous. In fact, the occurrence of a fairly abrupt change 
in one or all flow characteristics at a critical spacing is 
one of the fundamental observations of flow interference in 
cylinder arrays. 

For the tandem arrangement (one cylinder behind the other, 
see Figure 1), it has been shown that at relative spacings 
S/D<3.5 there is a strong mutual interference between the two 
cylinders. This critical distance increases with the bluffness 
of the bodies. For two plates in tandem, the critical spacings 
is about 10 times the plate width / 6 /. In general, the 
tandem arrangement has a strong effect on the drag and is sen- 


Sitive to spacing. The upstream cylinder takes the brunt of 


Low t 
aeeete Sp 
So ©) : 


—s-— 
Fig. 1 Tandem arrangement of two cylinders 
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the burden and the total drag for the group is smaller than 
the sum of the drag forces acting on each cylinder in isola- 
tion in a tunnel with the same blockage. The drag coefficient 
for the downstream cylinders is reduced partly by shielding 
and partly by the occurrence of earlier transition in the 


boundary layers due to the turbulence in the approaching flow. 


FLOW 


Fig. 2 Side-by-side arrangement of two cylinders 


The side-by-side arrangement (see Figure 2) exhibits a 
discontinuity in the flow and measured forces for spacing 
meeos smaller than a critical value (S/D<2.2 for two circu- 
lar cylinders and about 4 two plates). The observed dis- 
continuity (a switching phenomenon and bistable nature of 
the flow between plates) results in two values of the drag 
coefficient. This results from the mutual interference of 
the vortices on the adjacent sides of the vortex streets. 
The bistable nature of the flow may be avoided by increasing 
the size of one of the cylinders. The sum of the bistable 
high and low drag is often less than twice the drag of the 


Single cylinder. For two plates in side-by-side arrangement, 
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the said sum may be about 10 percent larger than twice the 
drag of the single plate. This increase strongly depends on 
the width of the flow field relative to the total width of 

the bodies, i.e., on the blockage ratio. For relatively 

large blockage, the flow is forced through the openings 
between the cylinders rather than around the total configura- 
tion. Consequently, the use of ordinary shielding and block- 
age factors for groups of cylinders is meaningless. Experi- 
ments must be conducted in channels or tunnels with very small 
blockage in order to obtain valid force-transfer coefficients. 
There is at present no means to separate the flow interference 
effects from the blockage effects for groups of cylinders. 

Additional work with tube arrays in steady flows has been 
reported by a number of people / 7 - 14 /. Ross / 15 / con- 
ducted large scale wave tank tests for the case of one cylin- 
der on each side of the test cylinder in the range of critical 
Reynolds numbers. His results appear to indicate that the 
wave force increases significantly only when the spacing 
between two cylinders is less than about one diameter. 

The foregoing studies, conducted mostly in steady flows, 
show that the sum of the individual drag forces for isolated 
cylinders is often greater than the total force acting on the 
array as long as the local boundaries do not constrain the 
flow. They further show that the results obtained from the 
tandem and side-by-side arrangements cannot be generalized to 
predict the combined interference resulting from the staggered 


arrangement. 
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Gibson and Wang / 16 / carried out two different experi- 
ments to determine the added mass of a series of tube bundles. 
The bundles consisted of tubes of uniform diameter d arranged 
either in a square configuration or a circular configuration. 
In the first series of experiments, they towed the model of a 
pile cluster under linear acceleration. In the second series, 
they have vibrated the model at its own natural frequency. For 
both cases, they have calculated the added mass through the 
use of the measured force and acceleration and plotted them as 
a function of the 'solidification ratio' defined by td. / 7D 
where D is the pitch diameter of the bundle. Their results 
have shown that the added mass increases sharply after the 
solidification ratio reaches the value of 0.4 to 0.5. Beyond 
this value, the volume enclosed rather than the volume dis- 
placed by the structure becomes important. This result is 
disputed lam] / and the results of both series of tests are 
no more applicable to separated wavy or oscillatary flows 
about tube bundles than those predicted from the potential 
theory with or without diffraction effects. 

Relatively few studies have been carried out with oscil- 
lating tube bundles /18, 19 7. Bushnell / 19 / determined 
the interference effects on the drag and transverse forces 
acting on a single member of a two-cylinder configuration 
through the use of a pulsating water tunnel. He did not 
evaluate the drag and inertia coefficients through the use 
of a suitable method, e.g. Fourier averaging. Instead, he 


picked out the maximum force values which occurred in each 
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half cycle and averaged them over ten consecutive values so as 
to obtain a mean maximum force for each flow direction. The 
results have shown that the presence of neighboring cylinders 
significantly affects the forces on an individual cylinder of 

an array and the interference effect increases with increasing 
relative flow displacement. The maximum drag force on shielded 
cylinders was reduced relative to an exposed cylinder by up to 
50 percent. Bushnell has suggested that a design using a high 
Reynolds number single-cylinder drag coefficient applied through- 
out the array would have an extra margin of safety against maxi- 
mum drag loading due to interference effects. The transverse 
force could be 3 to 4 times larger for interior array positions 
than that for a single cylinder. Thus, a cylinder array, such 
as a riser, supported at regular intervals may exhibit very 
complex dynamic behavior. Some of the members between supports 
may undergo violent transverse oscillations. 

It is evident from the foregoing that there is very little 
systematic data for tube bundles and two cylinder interference. 
In view of this fact is was decided to carry out experiments 
meen tCwo bundle configuration consisting of 12 and 9 outer 
tubes, respectively. Subsequent to this investigation a 
detailed study of the interference effects of two cylinders 


in harmonic flow was undertaken. 
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TIl. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A. OSCILLATING FLOW TUNNEL 

Experiments were carried out ina large U-shaped oscillating 
water tunnel (see Figure 3). It has been used extensively at 
Pataministitution over the past 6 years / 20 /. Only salient 
Peaeures, Most recent modifications, as well as the adaptation 
for this work, are briefly described here. 

The length of the tunnel has been increased from 30 feet to 
35 feet and its height from 16 feet to 24 feet. Previously, a 
butterfly-valve arrangement at the top of one of the legs was 
used to initiate the oscillations. In 1978, the tunnel was 
modified so that the oscillations could be generated and 
maintained indefinitely at the desired amplitude. For this 
purpose the output of a 2 HP fan was connected to the top of 
one of the legs of the tunnel with a large pipe. A small 
butterfly valve, placed in a special housing between the top 
of the tunnel and the supply line, oscillated harmonically 
at a frequency equal to the natural frequency of the oscilla- 
tions in the tunnel. The oscillation of the valve was per- 
fectly synchronized with that of the flow through the use of 
a feedback control system. The output of a pressure trans- 
ducer (sensing the instantaneous acceleration of the flow) 
was connected to an electronic speed-control unit coupled to 
a DC motor oscillating the valve plate. The circuit main- 
tained the period of oscillations of the valve within 0.001 


seconds. The fluid oscillated with a period of 6.000 seconds. 


ey 
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The amplitude of oscillations was varied by constricting or 
enlarging an orifice at the exit of the fan. The flow oscil- 
lated at a given amplitude as long as desired. 

Two identical force transducers, one at each end of the 
cylinder, were used to measure the instantaneous in-line and 
transverse forces. A special housing was built for each gage 
so that it can be mounted on the tunnel window and rotated to 
measure either the in-line or the transverse force alone. 

The calibration of the each gage was accomplished by 
hanging loads in the middle of the cylinder after setting both 
gages to sense only the transverse force (here in the vertical 
direction). Repeated calibrations have shown that (a) the 
gages were absolutely linear up to 500 pounds; (b) the gages 
yielded the same signal for loads applied either upward or 
downward; and that (c) the gages, together with the electronic 
system to which they are attached, were capable of sensing 
loads as small as 0.02 pounds. 

The displacement, velocity, and acceleration of the fluid 
in the tunnel were obtained through the use of several methods 
described in / 20 /. 

The data were recorded in analog form and then digitized 
and processed through the use of an HP9845 computer-digitizer 


system. 


fee LEST BODIES 
The tube bundles were made of aluminum rods and mounted 
between two circular stainless steel plates. Two configura- 


tions were chosen for test purposes. These configurations 
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roughly correspond to the most commonly used riser shapes in 
the off-shore industry. Figures 4a and 4b show the two riser 
configurations together with their characteristic dimensions. 

The two cylinder configuration was arranged in a similar 


manner (see Figure 5). The active cylinder was mounted on the 


Fig. 5 General arrangement of two cylinders 


Force transducer as described previously. The passive cylin- 
der spanned the entire test section and was rigidly mounted 
on the tunnel walls. Both cylinders were smooth and 3 inches 
in diameter. The angle a was varied from zero to 90 degrees 
in 30 degree increments. The N in cylinder spacing ND was 
varied from 1.5 to 3.5. The in-line and transverse forces 
acting on the active cylinder were measured and analyzed in 
terms of the appropriate lift coefficients and the drag and 
inertia coefficients through the use of the Fourier analysis. 
The details of the definition of the various force transfer 


coefficients are described in the following. 
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IV. FORCE TRANSFER COEFFICIENTS 


A. TUBE BUNDLES 
The in-line force, based on Morison's equation is written 


as 





Be —0.5/C ot Oy | Coswt| Coswt 
dm i. 
(i) 
3= (0) 5 Zs sea D<u wSinwt 
m im 
in which Un is the maximum velocity in a given cycle. 
Denoting 
D ==ID en 
a SL e al 
and inserting in equation (1), one has 
1m a 
=== =i -C4° (p=) [Cos wt | -Cos ut 2) 
ie5,D LU 2 
e°--m 
ne 
+ - Si 
TP a Sinwt 
(—) 
D 
e 
in which T = 27m. 
For an oscillating flow represented by U = -U. coswt the 


averages of C,. and Cc, are given by / 20 / 


d 
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3 e 277 
cee — - (—) ff ae. Oe op 
d 4 aa FAD UI 
e-m 
and 
2U_, At oe PF Sinwt 
Cc = a i) d ut (4) 
= t D O 2s De 
2 em 


Pine which ae represents the measured force. 


The lift coefficient is defined as 


| in which F, represents the amplitude of the transverse force 


L 
and L, the length of the cylinder. 
The vortex shedding frequency for the entire bundle was 


expressed in terms of the Strouhal number defined by 





in which f. represents the frequency of vortex shedding; T, 
the period of flow oscillation; and K, the Keulegan-Carpenter 


number defined by 
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B. TWO CYLINDER ARRANGEMENT 

The in-line force acting on the active cylinder is not in- 
dependent of the direction of the fluid motion because of the 
interference effect of the passive cylinder. In the interval 
t/2<ut<37/2 the direction of flow is opposite to that in the 
interval 3m/2<ut<57r/2. With this fact in mind Morison's equa- 


tion was written as 


=m ° C,Sinwt - C,/Cos ut|Cosut (6) 


and - and Cy were determined from 


4-U_-T es: 
a a Ru f 31/2 Ff oe ACS) (7) 
TF 1/2 fu, LD 
a =-5 f 31/2 a= di( wt) (8) 
7 2 ?fU_ LD 
m 
_ 8 Su/2 F.Sinut (9) 
2 =a 3 Lan/2 pen ae e d (wt) 
‘eID. Face G10, 
m 
and 
C oS f 5717/2 F-Cosuwt Se (10) 


Bin y 2 fU_, LD 
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in which Cal and Cal represent the inertia and drag coefficients 


mametne interval 1/2<wt<31/2; and C2 and C427 


coefficients for the interval 37n/2<ut<5r/2. In general these 


the corresponding 


two sets of in-line force coefficients are not identical accept 
when. a=90°. 

The evaluation of the transverse-force coefficients re- 
quires a distinction not only in the direction of flow but 
also in the direction of the transverse force. As the flow 
proceeds in a given direction the lift force acting on the 
active cylinder does not remain symmetrical primarily because 
ef the proximity of the passive cylinder, situated asymmetri- 
cally relative to the direction of the flow. It is also re- 
cognized that the random nature of the vortex shedding does 
not allow a perfectly symmetric lift force even in the absence 
of any interference effects. Consequently, one must define 
Bew@melitt coefficients. In the interval 7/2<uk<31/2, one 


has (see Fig. 6a) 


Lio .5/u* pr 
m 


and in the interval 31/2<wt<5T1/2, (see Fig. 6b) 
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1S) 
f3 ~ 0.5°U5 DL 
and 
Bead 
“L4 7 0.5/U~ DL 


where F through F 


bl are shown in Figures 6a and 6b. 
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The Strouhal number was defined in amanner similar to that 


previously described and evaluated through the use of 


ri 
Se. = 7 
m 
mewiech i = l, .. ee inedtcatdnme scene lift force directions 


Shown in Figures 6a and 6b. 


C. GOVERNING PARAMETERS 

It has been shown previously / 20 7 that the Fourier- 
Seemweaged lift, drag, inertia coefficients for a smooth cylin- 
der immersed in a harmonically oscillating flow depend on 
the Reynolds number and the Keulegan-Carpenter number. [In 


other words 


In terms of a frequency parameter 8 defined by 
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one has 





For the two cylinder arrangement one has to introduce two 
additional parameters namely, the angle a of the line joining 
the center of the two cylinders and the relative distance 
between the cylinders. Thus, all the force-transfer coeffi- 
cients cited previously depend on the Keulegan-Carpenter 
number, Reynolds number (or the frequency parameter), the 
angle a, and the relative distance factor N (see Figure 5). 


Thus one has 


C, (a coefficient) = EF. (K,8,a,N) 
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Vo beeecUsoLON OF RESULTS 


eee tbe BUNDLES 

The drag and inertia coefficients for the two riser con- 
figurations are shown in Figures 7a and 7b as a function of 
the Keulegan-Carpenter number for constant values of 8. The 
drag coefficient decreases gradually with increasing K. The 
inertia coefficient increases with increasing K and reaches 
values as large as 6. The use of the potential flow theory 
for unseparated flow yielded / 21 / an average inertia 
coefficient C* eel LOm ene EwOuCOntigubat ions . ome was 


defined by 


2 


2 
* = 
cS a (C.D) /2D. 


The comparison of co with that obtained experimentally has 
Shown that some fluid mass is entrapped within the array and 
that neither the potential theory nor the diffraction analy- 
Sis can adequately describe the behavior of the complex 
separated flow through the bundle. 

A comparison of the Figures 7a and 7b shows that the drag 
and inertia coefficents as defined here are nearly indepen- 
dent of the particular arrangement and the number of tubes 


comprising the bundle. 
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Experiments with different tube arrangements [EP 7 have 
shown that the force coefficients are independent of the 
Reynolds number within the range of Reynolds and Keulegan- 
Carpenter numbers reported herein. 


The reason for the dependence of C, and Cn On Kes. Enoughe 


d 
to be the dependence of the interaction of the wakes of the 
outer and inner pipes. The vortices in the wake of a given 
cylinder loose about 70 percent of their strength within 10 
cylinder diameters. Thus, for small values of K, the vortices 
generated by a small tube at the center front of the bundle 
arrive at the central tube as weak vortices. Consequently, 
each tube behaves more or less as if it were independently 
subjected to a turbulent harmonic flow. As K increases, not 
only the turbulence level but also the interaction between 
the wakes of the various cylinders increases. There is a 
certain amplitude of oscillation beyond which neither the 
interaction of the wakes nor the increase of the turbulence 
Bevel affects the overall force acting on the bundle. A com- 
parison of the total drag force acting on the bundle with the 
sum of the drag forces acting on each cylinder in isolation 
in harmonic flow (at the corresponding K and Re values 
appropriate to each tube) shows that the former is about 10 
percent smaller. 

The lift coefficients, as given by Eqn. (5), for the two 
riser configurations are shown in Figures 7c and 7d. Apparently, 
both configurations yield nearly identical lift coefficients. 


As will be seen later the rapid decrease of the lift coefficient 
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with increasing K is in conformity with that for a single 
cylinder. As the Keulegan-Carpenter number increases, each 
member of the bundle sheds an increasing number of vortices 
with different strengths and frequencies. 

Furthermore the shielding or proximity effect of the tubes 
results ina highly confused vortex system. This in turn 
decreases the coherence of the regular vortex shedding, re- 
sulting in very low lift coefficients. Evidently, a greater 
understanding of the lift force requires the determination of 
the lift force acting on each member individually and the 
consequence of tube proximity. 
peeeemeaG AND INERTIA COEFFICIENTS FOR THE TWO CYLINDER 

ARRANGEMENT 

The drag and inertia coefficients for N=1.5, N=2, N=2.5, 
and N=3.5 (only for a=0) and for a=0°, a=30°, a=60°, and a=90° 
are shown in Figures 8 through 20. Each figure represents 
either the drag or the inertia coefficient for two different 
flow directions. For example, Figure 8a shows the drag 
Goefficients C,, and C,, for a=0° and N=1.5. The coeffi- 
cients Cal is for the flow directed from passive towards 
active cylinder. The coefficient Ca fsmionetie flow din— 
ected from the active towards the passive cylinder. [In 
addition, included in each figure is the drag or inertia 
Coefficient for the single cylinder of the same size. 


Figure 8b shows the inertia coefficients Ch and C2 ior 


Jl 


are related to 


O 


@=0 and N=1.5. The coefficients os and C, 


a 2 


the flow direction in the same manner as the coefficients 


Cal and Cao- 


Sab 





Figures 8a, 9a, 10a, and lla show that for the tandem 


arrangement corresponding to a=0° the Grageecoerrt relent Ca 


is always larger than C The difference between the two 


ool 


drag coefficients decreases with increasing N (larger spacing). 
This finding is in conformity with the measurements reported 
by Zdravkovich / 5 / for the tandem arrangement in steady 

flow. For relatively small spacings (smaller N) the wake of 
the passive cylinder is highly turbulent. This causes earlier 


separation on the active cylinder and hence in smaller Cal 


values. One may also express the same fact by stating that 
the active cylinder is shielded by the passive cylinder when 
the flow sees the passive cylinder first. As the spacing 


between the two cylinders increases, and Ca gradually 


eat 


approach each other and single cylinder value. Even for 


M=3.5, C is slightly larger than C because of the afore- 


o dl 


mentioned effect of turbulence in the wake of passive cylinder. 


The inertia coefficients om and om for the tandem arrange- 


i 2 
ment are shown in Figures 8b, 9b, 10b, and 1llb. Unlike the 


drag coefficient, the inertia coefficient oe is always 


iz 


smaller than ce This is not entirely unexpected. All drag 


L: 
and inertia evaluations for single cylinders in waves and 
harmonic flows have shown that / 20 / the drag and inertia 
coefficients follow opposite trends. In other words, when 
the drag coefficient increases the inertia coefficient de- 
creases and vice versa. The reason for this is embedded in 


the formulation of the Morrison's equation (See Eq. (1)). When 


the phase angle between the velocity and the maximum force 
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meena the total force is primarily due to the velocity- 
square dependent drag and the contribution of the inertial 
force is negligible. In fact such flows are termed drag- 
dominated flows. When the phase angle between the velocity 
and maximum force is large, the total force is primarily due 
to the acceleration-dependent inertia. Such flows are termed 
Memama-dominated flows. Thus, it is apparent that the drag 
and inertia coefficients should follow opposite trends. 

As the relative spacing increases (See Fig. llb) Cal and 
C2 approach each other. However, the difference between two 
coefficients remains relatively larger than that between C 


dl 


and Cao° The experimental fact is that the spacing between 
the two cylinders in tandem arrangements must be considerably 
larger cere 3.5 diameters for the inertia coefficients to 
become independent of the proximity effects. It is also clear 
that structures in close proximity can cause relatively more 
Significant changes in the inertia coefficient than in the 
drag coefficient. 

Figures 8b, 9b, 10b, and llb also show that the inertia 
coefficient for single cylinder (dashed lines) is nearly the 
same as Cn2° In other words, when the flow sees the active 
cylinder first the proximity effect on the inertia coefficient 
C2 is negligible. It is also apparent that when the flow 
sees the passive cylinder first the inertia coefficient C41 
1s considerably increased, The important conclusion to be 


reached on the basis of these results is that the inertia 


coefficient of a shielded cylinder may increase considerably 
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and the drag coefficient may decrease significantly. 

The other extreme of the two cylinder arrangement is the 
Side-by-side arrangement. The discussion of this case will 
be taken up next since the force coefficients for all other 
intermediate values of a follow trends intermediate to the 
tandem and side-by-side arrangements. 

Figures 18a, 19a, and 20a show Cal and Ca fOrmN=1.5, N=2, 
N=2.5 for the side-by-side arrangement. Several facts are 
immediately apparent from these figures. For spacings as small 
as 1.5 diameter the drag coefficients are not greatly in- 
ane ae 


“a1 a2 
are very close to that for the single cylinder. The slight 


cluenced by the proximity effects. Secondly, 


difference in the lower range of Keulegan-Carpenter numbers 
May be attributed to the increased blockage in the tunnel and 
experimental uncertainty which was found to be about 5 per- 
cent. Experiments / 5 / with steady flows about the same 

type of two-cylinder arrangement have also shown that the 
proximity effects are negligible for spacings larger than 
about 1.5 diameter. 

shown in Figures 18b, 


mine inertia coefficients Ch and Ch 


i. Z 
19b, 20b exhibit a somewhat different behavior. Even though 
they are nearly identical for any given spacing they are con- 
Siderably larger than that for the single cylinder. Only for 
N=2.5 that they approach the single cylinder values. 

The two other arrangements intermediate to those discussed 


above behave as would be expected. Both the drag and inertia 


coefficients shown in Figures 12 through 17 vary between the 
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cases for a=0° and a=90° and both gradually approach the 
corresponding single cylinder values. The conclusions re- 
garding the effect of shielding remain valid for the inter- 


mediate arrangements. 


C. LIFT COEFFICIENTS FOR THE TWO CYLINDER ARRANGEMENT 

As in the case of drag and inertia coefficients, first the 
tandem arrangement, then the side-by-side arrangement, and 
Finally, the intermediate arrangements will be discussed. 

For the tandem arrangement C_,=C and C =Cr 4 (See Fags. 


el erie 2 ibys 


6a and 6b). The coefficients C Om ae are shown in Figures 


Ll L2 


21-24 and the coefficients C.3 or Cha are shown in Figures 
25-28. Also shown in each figure is the lift coefficient for 
a single cylinder of the same size. 

The data showed that the lift coefficient is strongly 
affected by the cylinder proximity, irrespective of the 
direction of the flow. For Keulegan-Carpenter numbers larger 
than about 20 the lift coefficients Cri and C.3 remain prac- 
tically constant. Only for very small Keulegan-Carpenter 
numbers that Cra and Cr 3 Inerease maoldiky, in Comfeormity with 
mite £LOr a single cylinder. The reason for this is that for 
Keulegan-Carpenter numbers near 10, at most two vortices are 
shed during a given half-cycle. These vortices do not travel 
1.6 diameters in one direction as seen from the simple relation 
that K=2TA/D where A is the amplitude of the flow oscillation. 
Consequently, the vortices shed from the active cylinder are 
only partially blocked by the passive cylinder for K smaller 
than about 10. Nevertheless, the effect of cylinder proximity 


is ever present as seen from a comparison of the data with 
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that for the single cylinder. As the spacing increases, the 
ieeeeecOecrtricients gradually approach that of the single cylin- 


der as seen in the Figures 24 and 28. 


In the case of side-by-side arrangement C11 =Cr3 and 
C27 era: The coefficients Chi OF Cr3 are shown in Figures 29- 
31 and the coefficients C or Char in Figures 32-34. The 


2 


lift coefficients for this particular arrangement have been 
measured for only N=1.5, N=2, and N=2.5. 

As the flow oscillates about the cylinders arranged in 
Side-by-side position, the blockage introduced by the proximity 
of the cylinders forces the fluid to flow above the top and 
below the bottom cylinder. In other words, the velocity on 
the extreme sides of the passive and active cylinder increases 
while the velocity between the two cylinders decreases. This 
in turn increases the pressure on the sides of the cylinders 
nearest each other and decreases the pressure on the opposite 


Sides. Consequently, C Cpe 16 Should be smaller than C 


ity i563 Hie 
or Cra: A comparison of the Figures 29-31 with Figures 
32-34 shows that Cri Or C13 1s in fact, always smaller than 
Ch2 Or Cha: As the relative spacing increases all four 


coefficients gradually approach that for a single cylinder as 
seen from Figures 31 and 34. These results are in conformity 
with those obtained with steady flow about a side-by-side 
arrangement of circular cylinders. / 5 / 

For arbitrary arrangements one must deal with four dif- 


ferent lift coefficients, namely, and Cc. 


Si! Sap? pee A. 


The coefficient Cri is shown in Figures 35-37 for a=30° and 
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for N=1.5, N=2, and N=2.5. Similarly; C 5 is shown Figures 


L 
38-40; C3) in Figures 41-43; and finally, Cha? in Figures 
44-46. The corresponding lift coefficients for a=60° arrange- 


ment is shown in Figures 47-58. 
A comparison of Figures 35-37 with Figures 38-40 shows 
iia t Ce 


is always larger than C The reason for this is 


1 poe 
that the vortices shed from the passive cylinder pass over 
the top of the active cylinder and thereby reduce the pressure 
at the top of the active cylinder. The coefficients Cr and 
Ch gradually approach that for a single cylinder as the 
Spacing increases. 

When the flow first sees the active cylinder, C3 remains 


relatively larger than C. because of the presence of the 


4 
vortices near the top of the active cylinder which were shed 
during the previous half cycle. As the spacing between the 
and C 


two cylinders increases both Cr rapidly approach the 


S L4 
lift coefficient for a single cylinder. 
For a=60° arrangement all four lift coefficients are 
practically unaffected by the cylinder proximity even for a 
spacing as small as 1.5 diameters and they are nearly equal 
to that of an isolated cylinder. Evidently, the effect of 
eyeunder-cylinder proximity on the lift coefficient for a 
harmonically oscillating flow is confined to relatively small 
spacings and angular positions. 

The frequency of the lift force oscillations has been 


expressed in terms of the Strouhal number and evaluated for 


all cylinder arrangements. The results have shown that the 


ay 








Strouhal number remains constant at an average value of about 
0.2. The only exception to that was the case corresponding 
to a=30° and N=1.5. 

Figures 59 and 60 show the Strouhal number as a function 
of K for a=0° and N=1.5. As N increased the Strouhal numbers 
for the two flow directions rapidly approached each other and 
a mean value of 0.2. 

Figures 61 and 62 show the Strouhal number for a=90° and 
N=2. Once again the lift coefficients for the two directions 
of flow are identical and gradually approach a value of about 
ee eeeeenm the case of a=30°, the wakes interfere in an asymme- 
tric manner and the Strouhal number shows a strong dependence 
on the Keulegan-Carpenter number. Figures 63-68 show the 
Strouhal number as a function of K for various values of N 
for a=30°. Apparently, the Strouhal number decreases from 
values as large as 0.4 to about 0.15 as K increases. For 
spacings larger than about two diameters Strouhal number 


decreases to about 0.2. 


oo 





Vinee cONCEUSIONS 


An extensive investigation of the forces acting on tube 
bundles and arrangements of two cylinders has warranted the 
following conclusions: 

1. Experiments with two particular riser configurations 
have shown that the inertia coefficient is considerably 
larger than that predicted by the potential theory 
and shows that some fluid mass is entrapped within the 
bundle as a consequence of the solidification of the 
mube Configuration. 

2. For tandem arrangement of two cylinders the drag and 
inertia coefficients depend on both the Keulegan- 
Carpenter number and the relative spacing of the 
cylinders. As the amplitude of flow oscillation be- 
comes comparable or smaller than the gap between the 
two cylinders, the drag and inertia coefficients 
gradually approach those corresponding to an isolated 
cylinder, i.e., the wakes of the two cylinders do not 
interfere with each other. 

3. The drag and inertia coefficients for the side-by-side 
arrangement exhibit a similar behavior. For a relative 
Spacing larger than about 2.5 diameters, the cylinders 
behave as if they were independent. 

4. The quantification of the lift force for an arbitrary 


arrangement requires four different lift coefficients, 


ce 





depending on the direction of the flow and the direc- 
tion of the force. These coefficients have been eval- 
uated and it has been shown that the effect of cylinder- 
cylinder proximity is confined to relatively small 
Spacings and angular positions. 

The Strouhal number for practically all two-cylinder 
arrangements remained nearly constant at about 0.2. 
Only in the case of 30° orientation did the SEONG 
asymmetric interference of the cylinder wakes result 
in a Keulegan-Carpenter-number dependent Strouhal 
number, ranging from 0.4 to 0.2. However, as the 
spacing increased the Strouhal number again approached 


a mean value of 0.2. 
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